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Tuning hole mobility in InP nanowires 
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Transport properties of holes in InP nanowires were calculated considering electron-phonon interaction via 
deformation potentials, the effect of temperature and strain fields. Using molecular dynamics, we simulate 
nanowire structures, LO-phonon energy renormalization and lifetime. The valence band ground state changes 
between light- and heavy-hole character, as the strain fields and the nanowire size are changed. Drastic changes 
in the mobility arise with the onset of resonance between the LO-phonons and the separation between valence 
subbands. 



Semiconductor nanowires (NWs) are increasingly used 
in a wide range of devices. They appear as building 
blocks of nanocircuits 1 and can be applied in electri- 
cally driven lasing 2 , which can be used in telecommunica- 
tions and information storage for medical diagnostics and 
therapeutics 3 . Improvements in NW synthesis, including 
chemical technics, allow thorough control of their shape, 
size and composition 4-8 along with detailed microscopic 
characterization of built-in strain fields. 8,9 As the con- 
ductivity is mostly defined by the carrier-phonon inter- 
action and phonon-lifetime, tuning of the NW structural 
properties could result in the possibility of also finding 
optimal conditions for carrier transport. 

Considerable efforts have been devoted to the descrip- 
tion of carriers in the conduction band of NWs, 10-15 while 
similar endeavors are not so common for holes in the va- 
lence band. 16 As the mobility is inversely proportional 
to the carrier effective mass, one may naturally expect 
that considering carriers in the valence band may result 
in a drop in mobility when compared to the light elec- 
trons in the conduction band. This could certainly be 
the case for heavy-hole (hh) transport; however, light- 
holes (lh) under certain conditions may be promoted to 
be the top valence band by tuning structural parame- 
ters of NWs. 17 This atypical circumstance is the result 
of confinement effects and hh-lh mixing, affected as well 
by strain and surface asymmetry fields. 8,17 As we will 
show here, this results in significant mobility enhance- 
ment for lh in suitable NWs. We can also take advan- 
tage of valence band mass anisotropy to attain resonant 
conditions that allow sharp variations of the hole mobil- 
ity with external parameters, especially when the leading 
scattering process involve longitudinal optical phonons 
(LO-phonons) through the deformation potential. 18 Ad- 
ditional hole-phonon interactions, 19 such as deformation 
potential and piezoelectric coupling to acoustic phonons 
and polar coupling to optical phonons, 20 have weaker ef- 
fects and will not be consider here. 19,21 In order to pro- 
vide realistic estimates of the expected mobility changes 
in the NWs of interest, we consider the effects of dimen- 
sionality reduction on the LO-phonon dispersion and life- 
time, using molecular dynamics simulations for different 
NWs size and at various temperatures. 



We consider different NW cross sections and shapes, 
while temperature effects are included in the mobility 
through the phonon occupation and strain effects in a 
multiband Luttinger Hamiltonian. We find that mobil- 
ity changes in a non-monotonic fashion according to NW 
width, strain fields, and temperature. In particular, we 
show that for certain NW widths, one finds resonant be- 
havior that greatly suppresses the hole mobility and is 
strongly affected by temperature and strain. The inter- 
action potentials used in our molecular dynamics (MD) 
simulations consist of two- and three-body interaction 
terms, as described by Branicio et a/. 23-27 . The param- 
eters of the interatomic potential are determined using 
the cohesive energy, density, bulk modulus and clastic 
constant C\\ of the material as described before, 27 with 
some slight adjustments. 28 

This potential provides excellent estimates for melting 
temperature, structural phase transformation induced by 
pressure, and specific heat 24 and describes well the vibra- 
tional density of states of the material. We obtain the 
phonon density of states in InP NWs considering differ- 
ent temperatures. The NWs were created cutting a block 
of a perfect crystal with the z-axis along the [001] direc- 
tion with periodic boundary conditions in the z-direction. 
The x- and y-directions were surrounded by a vacuum 
region. The system consisted typically of nine unit cells 
along x- and y-directions and forty unit cells along the z- 
direction (53Ax53Ax234.5A); the total number of atoms 
is 25,920 (12,960 In + 12,960 P) (Fig. 1(d)). The NW 
is allowed to relax during a long simulation run (25,000 
time steps, one time step=1.5fs) at each temperature. 
After this relaxation time, a few surface defects can be 
observed. 

The velocity- velocity auto correlation function, Z (t) = 
^(^'^.(o)) ' wnere Uj (t) is the velocity of particle i at time 
t and the brackets are averages over ensembles and par- 
ticles. The vibrational phonon density of states (VDOS) 
is determined through the Fourier transform 22 
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where the subindex I is the atom In or P. 
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FIG. 1. VDOS for NW and bulk, (a) Total VDOS; rigid ion 
model (RIM) VDOS for bulk from data. 30 (b) and (c) Partial 
VDOS for bulk and NW at 300if : (b) indium contribution 
and (c) phosphorus contribution, (e) and (f) VDOS for bulk 
and NW at 10 K and 300 K. (d) Simulated InP NW structure 
by molecular dynamics at T = 300if. Green (grey) dots 
represent Indium (Phosphorous) atoms. 



Fig. 1 (a) compares the VDOS computed from MD cal- 
culations for a bulk system (solid curve) to the VDOS 
extracted from the rigid ion model (dashed curve) based 
on experimental results 29,30 . The MD results reproduce 
very well the main characteristics of the experimental, 
predicting the existence of the transversal optical (TO, 
~ 43meV) and longitudinal optical (LO, <~ 38meV) 
modes, and a gap between 22meV and 36meV (acoustic 
modes below 22meV not shown). 

The effect of the surfaces on the VDOS of the NW is 
also shown in Fig. l(a)-(c). The main NW character- 
istics in the VDOS resemble the bulk results; however, 
some differences can be observed. The NW surfaces ap- 
preciably increases the amount of modes in the gap re- 
gion, between 22meV and 36meV. Although the TO 
and LO modes are dominant, the NW surface plays an 
important role by inducing surface modes that appear in 
the gap region. We notice that the P atoms make the 
main contribution to the optical modes in the gap have 
a larger VDOS, and an overall large mode width for the 
NW. The optical modes arc characterized by the relative 
displacement between ions and the lighter atoms usually 
dominate such modes, as shown. Moreover, we observe a 
slight blue shift of the LO mode (~ lmeV) with respect 
to the bulk. 

The effect of temperature on the VDOS for both NW 
and bulk structure is shown in Figs. 1(e) and (f). We 
see a general broadening of the modes and a shift in the 
peak position to lower frequencies with increasing tem- 
perature. Although the gaps region keep a similar profile, 
the VDOS in the gaps increase with temperature. We 
notice also that the temperature affects the contrast be- 
tween LO and TO modes, especially in the NW, Fig. 1(f), 



by increasing the LO amplitude. 

Notice that the mode broadening is slightly weaker for 
the NW than for the bulk. Based on the MD results, we 
verify that the TO and LO modes are the most important 
vibrational modes. Also, the NW exibits a shift in the 
LO mode with temperature, which is an important fact to 
be considered when calculating the contribution of hole- 
phonon scattering to the mobility of NWs. 

The valence band Hamiltonian for the NW can de- 
scribe confinement effects, mass anisotropy, and strain 
fields within the same framework, as one writes, 31 
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for the heavy- and light-holes, where j a (a — 1,2,3) 
are the Luttinger parameters, {A,B} = ^(AB + BA), 

and k± = k x ± ik y . Notice that the subband with hh 
character along the wire has a low effective mass in the 
transverse direction w (71 +72)^, while the lh subband 
has a large transverse mass sw (71 — 72) , the differ- 
ent transverse masses result in the possible inversion of 
the lh and hh subband ordering, due to the NW con- 
finement effects. Strain effects lead to modulation of the 
valence subbands, 8 introducing a subband displacement 

given by 32 AHhh = -P + Q and AHih = ~P + Q + ^, 
where P = 2(a v + a c )(^f^)e ll , Q = -b{^±^)e\\ 
and A so = Q.lQ%eV is the spin-orbit split-off energy. 33 
The hole wave function in the NW has the form |\ffj) = 
\ij}i)\J,mj), where \ipi) is the envelope function, which 
depends on the cross section of the NW, and \J, rrij) is 
the total angular momentum eigenstate, [3/2, ±3/2) for 
pure hh character, and |3/2, ±1/2) for the lh. 

The holc-phonon interaction Hamiltonian is given by 
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with M q = (q • e q )(h/2pLijqV) 2 , where q is the phoiion 
wave vector for polarization vector e q , p is the mass den- 
sity, V is the system volume, and Uh~i{q) is given in 
terms of the deformation potential for holes. 21,36 Consid- 
ering long wavelength processes, we have E//j_j(q) oc u, 
where u is the relative displacement between atoms in- 
side the primitive unit cell. By symmetry, states with 
hh character couple with those of lh character along the 
direction [001] 21,36 . As this coincides with the wire axis, 
along which the carrier transport takes place, we find 
(hh^\U h -i(z)\lh^) = (IfFph-iWhh*) - where 
do is the deformation constant and clq the lattice param- 
eter. 38 

The mobility is given by /1 = t, in terms of the 

hole-phonon scattering time, r _1 = ^2 a S(k,q), and the 
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transition rate 

S(k, k>) = — [\(Ul_ p )\ 2 S{E f (k') - Ei(k) - fiw q )+ 

\{H e h _ p )\ 2 5(E f (k')-E i (k) + hu (l )] (5) 

where k and fc' refers to the initial and final states, and 
"H-h^p and T-L^-p refer to the phonon absorption and emis- 
sion processes in Eq. (4). The phonon density is assumed 
to be given by a Lorentzian centered at lulo with width 
7. Both of these values shift with temperature, as dis- 
cussed in the previous section. We may now analyze the 
effects of strain and temperature on the hole mobility. 
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FIG. 3. Hole mobility versus NW width W for states with 
k z — 0. (a) System at T= 70K without strain, (b) System at 
T= 300K without strain, (c) System at T= 70K with strain, 
(d) System at T= 300K with strain. 



FIG. 2. Valence band ground states for NW of width W. (a) 
When the lh occupies the ground state, (b) When the hh 
occupies the ground state. 

To characterize the initial and final states involved 
in the scattering processes that affect the mobility, we 
show in Fig. 2 the relevant valence band structure for 
two different cases. For thin NWs, with or without 
strain, the finite NW width leads to a picture similar to 
Fig. 2(a), where the lh subband is promoted to the top 
given its higher transverse effective mass, as discussed 
before. Thus, under such conditions, a hh can be scat- 
tered to the subband with lh character through phonon 
emission (process Ei), and at T > the lh can be ex- 
cited to the hh subband via phonon absorption (process 
Ai). In the presence of lateral compressive strain, the 
subbands may switch their relative positions with the hh 
assuming the top at large NW width. Then, a lh might 
be scattered via phonon emission (process E2) while a 
hh can be affected by phonon absorption at T > (pro- 
cess A2). Notice that by changing the wire radius, one 
can reach a resonant condition (A£„j, = Tiujlo)- On the 
other hand, with strain, depending on the value of the 
NW width, the ground state can have a character lh (thin 
NW) or hh (thick). This behavior is similar for all NW 
cross sections, as it reflect the transversal quantization. 

The relative position of the valence subbands is ex- 
tremely important for the carrier transport in NWs. 
Given the mobility dependence on the longitudinal effec- 
tive mass, which modulates the hole-phonon interaction, 
valence subband shifts may produce sharp fluctuations 
of the mobility as temperature or structural parameters 
change. Fig. 3 shows the mobility for different strain 
and temperature values as function of the NW width. In 
Figs. 3(a)-(b), the mobility reflects a band configuration 
similar to the one depicted in Fig. 2(a). At T= 70 A', in 
Fig. 3(a), the increase effect of phonon absortion leads to 
the monotonic decrease of the lh mobility whit increasing 



NW width, while the hh displays a monotonic mobility 
increase, as the intersubband separation decreases with 
increasing NW width. Also, a sharp variation near the 
region where AE v i ~ Tt-^lo is seen, the resonant con- 
dition greatly enhances LO phonon emission by a hh in 
panel (a). At higher temperatures, Fig. 3(b), the reso- 
nant condition also affects the carriers in the lh subband, 
producing a sharp drop in mobility (from Ai processes). 

Given the band structure modulation with strain, the 
condition AE^bl ~ fu^LO can be attained twice by vary- 
ing the NW width (corresponding to the cases displayed 
in Figs. 2(a) and (b)). Thus, two resonant regions appear 
in Fig. 3(c) where phonons can be emitted by both the 
hh and lh subbands (Ei an E2 processes, respectively). 
At higher temperatures, the phonon absorption features 
appear as additional jumps in the mobility, shown in 
Fig. 3(d), processes Ai and A2. Notice that the lh and 
hh subband inversion with increasing NW width, in the 
presence of strain, is accompanied by crossing of the mo- 
bility curves, Figs. 3(c)-(d). Tuning the mobility of a 
hole system via in-situ changes of the NW width or strain 
fields, is not an easy task in experiments. As we will see 
below, however, one can achieve drastic in-situ mobil- 
ity changes for NWs close to the resonance condition by 
suitable changes in temperature. 

In Fig. 4(a) for a square NW width of W = 20A, no 
resonant signatures appear in the strain and temperature 
range analyzed, while for thicker NWs such conditions 
become evident for both kind of holes. This situation 
follows the trends described in Fig. 2, since for thinner 
NWs, the valence band configuration corresponds to that 
in Fig. 2(a), with AE v b > Hulo- In turn, by raising the 
temperature, the process of phonon absorption becomes 
more effective, reducing the lh mobility at higher temper- 
atures; as displayed in Figs. 4(b). The effect is present 
for all NWs, although with different features for various 
W values, as we now discuss. 
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FIG. 4. (a) and (b) hole mobility for two values of NW width 
for states with k z — 0: (a) as function of strain for T= 70K, 
and (b) as function of temperature for en = 0.9%. (c) Light 
hole mobility ratio, where Afiih = (fiih[70K] — {j,ih[300K]) , 
versus wire width for different values of strain at k z = 0. 



Fig. 4(c) shows the strongly non-monotonic mobil- 
ity variation with temperature for certain NW widths 
with energies close to the resonant condition, AE v b ~ 
hujLo- This figure shows the lh mobility ratio of high 
and low temperatures, A/xj/i//ij/,[70if], where A/i;/j = 
(M h [7QK} - m h [300K}), as function of the NW width for 
different values of strain. For lh-holes, the mobility at 
low temperatures is high (Fig. 4(b)) and according to the 
radius and strain one may observe a drastic drop in mo- 
bility at high temperature. For a system with no strain, 
for example, a large but monotonic drop in the mobil- 
ity is seen for NWs with large width (larger than lOOA). 
In real systems, however, free standing NWs grow with 
built-in strain, 8 which has a direct impact on the de- 
pendence of (i on the temperature and width. As the 
strain increases, the lh-hh subband reversal is possible 
as the NW width increases, resulting in highly sensitive 
mobility on temperature and/or width. For example, for 
£\\ = 0.9%, 8 and for width values close to 50 A mobility 
exhibit a sharp change in the mobility with temperature. 
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